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Abstract 



The finite-size scaling properties of the quantum Ising chain with different types of 
generalized defects are studied. These not only mean an alteration of the couphng 
constant as previously examined, but an additional arbitrary transformation in the 
algebra of observables at one site of the chain. One can distinguish between two 
classes of generalized defects: those which do not affect the finite-size integrability of 
the Ising chain, and on the other hand those that destroy this property. In this context, 
finite-size integrability is always understood as a synonym for the possibility to write 
the Hamiltonian of the finite chain as a bilinear expression in fermionic operators by 
means of a Jordan- Wigner transformation. Concerning the first type of defect, an 
exact solution for the scaling spectrum is obtained for the most universal defect that 
preserves the global Z2 symmetry of the chain. It is shown that in the continuum 
limit this yields the same result as for one properly chosen 'ordinary' defect, that is 
changing the coupling constant only, and thus the finite-size scaling spectra can be 
described by irreps of a shifted u(l) Kac-Moody algebra. The other type of defect is 
examined by means of numerical finite-size calculations. In contrast to the first case, 
these suggest a non-continuous dependence of the scaling dimensions on the defect 
parameters. A conjecture for the operator content involving only one primary field of 
a Virasoro algebra with central charge c = ^ is given. 



(This is a reprint version of a paper that appeared in Nucl. Phys. B 340 (1990) 633-658.) 



1 Introduction 



The one-dimensional quantum Ising chain, the Hamihonian hmit of the classical two-dimensional 
Ising model PP obtained via the transfer-matrix approach (see |2| E] and references contained 
therein) belongs to the most frequently investigated and best understood systems in statistical 
physics. The main cause of this is the fact that the Ising model (without an external magnetic 
field, however) admits an analytic solution ^ and on top of this undergoes a continuous phase 
transition at a finite temperature. It is known |S] that at criticality the Ising model describes the 
field theory of a free massless Majorana fermion, which is a conformal field theory with the central 
charge c of the Virasoro algebra (VA) being c = ^. 

After the importance of the conformal structure [HI [3 |H] of spin chains at the critical point 
had been recognized, the question arose how defects such that the system stays critical will affect 
these properties. Clearly one started the investigation of this problem with the simplest example, 
that is with the Ising quantum chain. The defects in this case were introduced by changing the 
coupling constant of one nearest-neighbour coupling corresponding to a change of the coupling 
constant on a half-infinite line in the two-dimensional classical model. The system in this case 
stays critical ^ and part of the conformal structure is maintained |10[I11[[T^I13| . 

It is the aim of this article to give a generalization of the case of one localized defect. This is 
done by performing an arbitrary transformation in the algebra of observables at one site of the 
chain. In other words, this means that one replaces the Pauli matrix in one coupling term 
of the iV-site Hamiltonian H of the (ferromagnetic) critical quantum Ising chain with periodic 
boundary conditions 

1 ^ 

^ = - 2 E + '^i ^J+i] ' ^N+i = , (1-1) 

by a general hermitian 2x2 matrix, therefore one has four real parameters ei, e^;, ej,, and that 
determine the generalized defect by 

c^w+i ==eil + ej,(Ti +ey(Tj'-|-ezCri , (1.2) 

where 1 denotes the 2x2 unit matrix and an 'ordinary' defect corresponds to ei—ey = ez~Q. 

The first results for the Ising quantum chain with 'ordinary' defects were obtained by detailed 
numerical calculation [TD|. but recently 113) the exact solution for an arbitrary number of 
isolated defects has been found. The algebraic structure of the spectrum for a fixed number of 
defects at commensurate distances has been examined in detail in ^2] (see also The results 

of these investigations for the case of one defect, where the spectrum generating algebra turns 
out to be a supersymmetric extension (the | = ^ higher symmetry algebra ^1]) of a so-called 
'shifted' m(1) Kac-Moody algebra (will be abbreviated by KMA in the sequel), are summarized in 
the subsequent section. Even more recently, another type of defect (a so-called 'extended' defect), 
where the coupling constant is altered on a whole fraction of the chain, has been studied [15) and 
a similar algebraic structure as above has been obtained. 

In contrast to the case of an 'ordinary' defect, a generalized defect as defined above might or 
might not destroy the finite-size integrability property of the chain and one expects a quite different 
behaviour for various types of generalized defects. Here and in what follows, finite-size integrability 
always means that one is able to write the Hamiltonian of the finite chain as a bilinear expression 
in fermionic operators via a Jordan- Wigner transformation [S] . A generic example consisting of a 
two-parameter family of defects for each of this two classes is examined by analytical respective 
numerical methods. In the case where the chain stays finite-size integrable the results of the exact 
solution resemble those of one 'ordinary' defect, which is contained therein for a special choice of 
the parameters. In particular, the scaling dimensions depend continuously on the parameters of 
the defects, whereas in the other case the numerical results favour a non-continuous behaviour. 
The operator content in the latter case is conjectured to be given by a single irrep of a VA with 
central charge c = ^ for fixed parameters of the defect. This article restricts to the quantum chain 
point of view, physical implications for the two-dimensional classical system — if any - are not 
discussed in this place. 
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The paper is organized as follows. In the subsequent section the previously examined boundary 
conditions and 'ordinary' defects together with the 'generalized' defects (|1.2ll are introduced more 
accurately. A short summary of the formerly obtained results is given. This includes a presentation 
of the character functions of unitary highest weight irreps of the VA with central charge c = ^ and 
of a shifted u{l) KMA. The generalized defects are then divided into two groups: one type that 
leaves the quantum chain finite-size integrable and the other that destroys this property. 

The third section contains an exact solution for the spectrum of the chain with an finite-size 
integrable type of defect that maintains the global Z2 symmetry of the Ising chain corresponding 
to 61=62 = in Eq. 1)1.2(1 . This is done employing the methods developed by Lieb, Schultz and 
Mattis (see also j^^^]) in a slightly generalized form. The operator content is given in terms 
of character functions of irreducible representations (irreps) of a shifted u{l) KMA and in fact 
turns out to have the same form as for one 'ordinary' defect )12) . whereas the behaviour on finite 
lattices appears to be different in general. 

In the forth part an example for a two-parameter class of finite-size non-integrable defects that 
corresponds to taking 61 =6j, = in Eq. H1.2|l is studied. This can be done through the application 
of numerical methods only. The effect of this type of defect on the finite-size scaling limit appears 
to differ qualitatively from the results of Sec. 3. The numerical data indicate that the dependence 
of the critical exponents on the parameters of the defect is not a continuous one. A conjecture for 
the operator content based on finite-size calculations for chains with a length of up to sixteen sites 
is presented. It can be expressed in terms of irreps of the VA with central charge c = i involving 
only one primary field in each case. 

The main results of this paper are summarized in the concluding Sec. 5. On top of this some 
questions that remained open are discussed in this place. The two Appendices consist of proofs 
concerning parts of the exact solution in Sec. 3 that had been omitted therein for the sake of the 
reader. 



2 'Ordinary' and 'generalized' defects 

In terms of the Pauli matrices cr^ , cr^ , and , the TV-site ferromagnetic quantum Ising Hamiltonian 
reads 

1 ^ 

^ = -2Eh^ + M-J+i] , (2.1) 
i=i 

where A plays the role of the inverse temperature and the meaning of has to be fixed 

by an appropriate boundary condition (BC). In this article only those BCs (or 'defects') will be 
considered that can be expressed in this way and where in addition o^^_^_l depends only on variables 
that live on the first site. Thus in Eq. H2.1|l only (at most) two adjacent sites are coupled directly. 
The Hamiltonian H2.1|l with a proper choice of the boundary term enjoys a global Z2-invariance. 
It commutes with the 'charge' operator Q defined by 

N 

where 1 denotes the identity operator. One can therefore divide the spectrum of H (|2.1|l into 
two sectors corresponding to an eigenvalue (—1)"^ of Q with ge{0, 1}. The system described 
by H (|2.1|) is the Hamiltonian limit of the classical two-dimensional Ising model that exhibits a 
continuous phase transition at A = 1 and the finite-size scaling limit of the spectrum is known for 
various different BCs, such as toroidal, free, fixed and mixed BCs ^1 |S1 El and for the case of 
an arbitrary number of isolated ('ordinary') defects ^2^]. Before starting with the subject of 
this paper and defining new possible choices of BCs, let us have a look to the formerly used BCs 
and the results that have been obtained for these. 

One already mentioned group of BCs is the so-called 'one defect' (10| (in the sequel the notion 
'ordinary' defect will be used to avoid confusion with the generalized defects) of strength Cx, where 



2 



one defines (J^^i to be (cf. Ea. (|1.2|) ') 

^Af+i — ■ (2.3) 

with e R, which in the two-dimensional model corresponds to a change of the coupling constant 
on a half- infinite line ||10j. For special values of the parameter e^, this includes most of the 
previously examined possible boundary terms. There are two distinguished choices for o^ff^i in 
(|2.1|) that correspond to periodic and antiperiodic BCs given by 6^ = 1 and e^ — — !, respectively. 
These are the sole BCs that allow the definition of a translation operator that commutes with H 
(|2.1|l and therefore are called toroidal BCs (see JHI)- Free BCs occur in this context as the case 
63; = 0. The also previously examined so-called fixed [Ifij BCs — which by performing a duality 
transformation can be shown to yield the same results as free BCs — as well as mixed BCs [T?i] 
cannot be written in this simple form, they are obtained by attaching constant terms to the two 
ends, respectively to one end (leaving the other end free), of the chain. 

All the above-mentioned BCs (or defects) have in common that the Hamiltonian H (I2.1|l 
commutes with the charge operator Q H2.2|l and beyond it can be written as a bilinear expression 
in fcrmionic operators. The last property ensures that the system stays finite-size integrable in 
each case. Therefore, the finite-size scaling limit spectra are known exactly for all these cases. In 
general, they can be expressed in terms of irreps of a 'shifted' u(l) KMA |14[ll9j that contains a 
VA with central extension c~l. For the special cases of free (as also for fixed and mixed) BCs 
it is also possible to use one VA with central charge c = ^ to describe the spectrum in the finite- 
size scaling limit (|H1E|) see also Sec. 4), whereas for toroidal boundary conditions one has two 
commuting VAs with c = ^ |S| due to translational invariance. 

To be more precise, let me reformulate the contents of the last paragraph in a more accurate 
way. As already mentioned, the operator content for free and toroidal BCs can be expressed |51ll7| 
in terms of unitary irreps of one respective two commuting VAs with c = ^ . For general central 
extension c, the Virasoro generators , toGZ, fulfill the commutation relations 

[Ljn, Ln] = (m - n)Lra+n + -^m{m'^ - l)S,n+ruO (2-4) 

with m,n E Z. The VA with central charge c=^ has only three unitary irreps [7] corresponding 
to a highest weight Ag{0, jg, ^}. The character functions xa{z) = tr(z^") in this case are given 
by EDI E] 

Xo(z) = ^ (z(12"^+")-z(12"^+7"+l)) •ny(z) = ^z4"'+".ny(z2) 

n6Z raeZ 



XJ-(^) 



nez 



^(12n^+2n) 



^(12n^ + 14n-|-4)^ 



Tlv{z) 



Hv{z) 



XI (z) = Z5 . ^ |^^(12n^+5„) _ ^(12n^ + lln+2)^ . = ,J . ^^4„=+3n . ^^,(^2) 

nez nez 

(2.5) 

Here, ny(z) denotes the function generating the number of partitions, which is given as an infinite 
product by the expression 

nW-) = ft (t^) • (2-6) 

m=l ^ ^ 

The partition functions T^{z, z) resp. ^q{z) for the finite-size scaling spectra of the Hamiltonian 
H H2.1|) in the charge sector q, qE{0, 1} are defined by |S] 



T^{z,z) = lim > 2 

'"'^ (2.7) 
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for toroidal BCs and |51[T71 



T„{z) = lim Vz^-'^W 

(2.8) 

for free BCs. Here, E^i,{P,N) denotes the k*^ eigenvalue of the Hamiltonian (|2.1|l with periodic 
(+) — resp. antiperiodic (— ) — BCs in the charge sector q and with the momentum (eigenvalue 
of the translation operator) P and the sum covers the whole spectrum of H 12.1|l . The analogue is 
true for the energies E^„.i^[N) for free BCs, except that in this case there is no translation operator 
commuting with H H2.1|l . The partition functions H2.7|l and 12.8|l are then given in terms of the 
character functions tabulated in Eq. (|2.5II by |SI 

= Xo(z)xo(z) + Xi(2)xi(^) 

T^iiz.z) = x^{z)xo{z) + Xo{z)x^{z) (2.9) 

Tq{z.z) = T+{z,z) = x-i^{z)x^{z) 

and by la ini 

Mz) = Xo{z) , Ti{z) - xi{z) . (2.10) 

For completeness, let me state the results for fixed and mixed BCs. The Hamiltonians for these 
BCs can be summarized in the form 

HihjN) = -\{jl H'^J+i + ^j] +^N + h-^l + lN- j (2.11) 

with /i, /at e {1, 0, —1} . Here, the different sectors of fixed BCs correspond to the choices /i =±1 
and /Ar = ±l, whereas mixed BCs are defined by /i = and /^r^il (or, equivalently /n — O and 
/i =±1). The partition functions for the finite-size scaling spectra in each case are defined relative 
to the ground-state of the chain with free BCs as in Eq. 12. 8|) . For fixed BCs 'S' , they are given 
by Xo(^) for equal fixed boundary spins (that is for _ff (1, 1) and _ff(— 1,— 1) ) and by Xi(^) for 
different fixed boundary spins (that is for H{1, —1) and H{1, —1) ) at the two ends of the chain, 
with these two BCs being related to the two charge sectors of free BCs by means of a duality 
transformation. In the case of mixed BCs [5j , one obtains XJ-(^) for the properly defined (that 
is, a non-universal surface term has to be taken into account) partition function of the finite-size 
scaling spectrum of H{0, ±1) H2.11|l . 

Now, let me come back to the general case of one 'ordinary' defect of strength e^. (|2.3|l . Instead 
of the Hamiltonian Hicx) H2.1|l with the defect (|2.3|) one considers the mixed-sector Hamiltonian 
H{e^) [121 

H{K) = H{e^)-Po + H{-e^)-P^ , (2.12) 
where Pq denotes the projection operator 

p, = \-{Q + {-m) (2.13) 

onto the eigenspace belonging to an eigenvalue (—1)^ of the charge operator Q (I2.2|l . {0, 1}. In 
this way one avoids the appearance of a non-local number operator |12| in the Hamiltonian. Of 
course, one can reconstruct the spectrum of H H2.1|l from the knowledge of the spectrum of H 
(|2.12f) . thus no information has been lost. The partition function 8{ex, z) of the finite-size scaling 
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spectra (relative to the ground-state energy for periodic BCs (cx — l), cf. for the chain with 
one defect of strength Cx is then defined by (cf. Eqs. H2.7|) and H2.8|) ') 



£(ex, z) — hm > 2 



(2.14) 

— N 

Ek{ex,N) = —■{Ek{ex,N) + Ai{ex)~E^{\,N)) , 

ZTT 

where Ek{ex,N) denotes the k*^ eigenvalue of the TV-site Hamiltonian H H2.12|l with the defect 
Cx (|2.3I) . The non- universal surface-energy term Ai{ex) that occurs in the expansion of the lowest 
eigenvalue Eo{ex) of the Hamiltonian 12.12|l with the defect ex |22l I23| 

- Eoicx) = Ao-N + Aiicx) + A2iex) ■ + o{N-^) (2.15) 

has been subtracted (the bulk term = § of course is not affected by a defect that alters only one 
term in the Hamiltonian (|2.12() ^. The dependence of this surface term on Cx is known analytically 
[T^ and it is given by 

1 + 1 
Ai{ex) = arctan(e:r) - - , (2.16) 

TTCx 2 

thence in particular ^i(O) = ^ — ^ and Ai{zLl) = 0, as necessary for toroidal BCs. Further- 
more, the values of A2{ex) in these cases are given by jS] A2(l) = 7r/12 , A2(— 1) = — 7r/6 , and 

A2(0) = 7r/48. 

The finite-size scaling limit spectra £{ex, z) (|2.14() can be expressed ^2 terms of character 
functions of unitary irreps of a shifted w(l) KMA. For this purpose define the tt(l) KMA by the 
commutation rules 

[TrmTn] = ITl ■ Sm+n,0 

(2.17) 

with 171,1% gZ , where the operators L,„ generate a VA (|2.4|l with the central charge c—1. The 
Lm result from the Tm by Sugawaras construction 51 through L„i — : Tm-rTr ■ , where the 

colons denote normal ordering. One then introduces a shift tp by means of the automorphism 

Tm I > Tm{f) = Tm + (p-Sm,0 

(2.18) 

Lm I > Lmif) = Lm+^-Tm+2^ ' ^mfi 

that leaves the commutation relations H2.17|l and ()2.4|) invariant. The unitary highest weight irreps 
of this algebra have the character functions 

Xt{z, y) = tr (z-^o^'^) • /o^'^)) = ziit+vf . n^(^) . yt+v ^ (2.19) 

indexed by a single real number <gR which is the eigenvalue of Tq. Here, Hy(z) is the function 
previously defined in Eq. H2.6|l . In the physical application one is interested in the eigenvalue 
of Lo{ip) in Eq. H2.19|l only, therefore one uses the character functions setting y=l and thereby 
leaving the actual value of t (there are in general two different ones) unspecified. 

The partition functions £{ex, z) l|2.14|l of the finite-size scaling spectra for the Hamiltonian H 
H2.12(l with one defect of strength Cx are then given by 



£{ex,z) = tr (^z^oM<=^'))j = ^ . n^(^) (3.20) 

nez 

with if{ex) determined by [111 I12| 

f{ex) = 7 - - arctan(ea:) . (2-21) 

4 TT 
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One recognizes that this is an infinite sum of characters of irreps of the tp - shifted u{l) KM A 
(|2.18|l . Equivalently one could interpret this as a single irrep of the shifted ^ — \ higher symmetry 
algebra P^ . 

After this short resume of known facts about the Ising quantum chain let me come back to the 
purpose of this article. As already defined in Eq. (|1.2|l of the preceding section, the notion of a 
'generalized' defect means that one introduces a boundary term by defining u%j^i in Eq. (|2.1|1 to 
be any expression subject to the sole requirement that it involves only variables living on the first 
site of the chain (and, of course, that the Hamiltonian (|2.1(l stays hermitian). The most universal 
term compatible with this condition will obviously be of the form H1.2|l . thus one has four real 
parameters that determine the defect. It is clear that all those defects with ei or Cz different from 
zero will destroy the global Z2 symmetry of the chain. Therefore, one can distinguish between 
those defects that preserve the symmetry and those, which break it (for more complicated models 
that the Ising model, there might be diff'erent degrees of symmetry breaking by defects that one 
would have to consider separately). More important for the analysis is the question, if the defect 
will influence the finite-size integrability property of the chain, thus this will be the classification 
of the 'generalized' defects that will be used in what follows. 

It is not the aim of this paper to exhaust all possible choices in Eq. (|1.2|) . but rather to 
restrict to the examination of two generic cases, each consisting of two-parameter families of 
generalized defects. In the subsequent section, the most general defect preserving the global Z2 
symmetry of the chain, that is ei = ez = in Eq. H1.2|l . will be studied. This always leads to an 
finite-size integrable Hamiltonian (however, it is not the most general type of defect doing this, 
another example for defects that share this property are for instance given by taking in Eq. (|1.2|l 
ea; = ej, = e2 = and ei arbitrary, which for ei^^O break the Z2 symmetry). A completely different 
type of defect is considered in Sec. 4. This is given by letting in Eq. (|1.2|l ei and Cy be equal 
to zero and the other two coefficients unspecified. The defects of this type have the property to 
destroy both the global symmetry and the finite-size integrability of the chain as long as does 
not vanish. 



3 The exact solution for one finite-size integrable general- 
ized defect 

Consider the two-parameter family of Hamiltonians 

ff(a,0) = -\ [a| + « + ^-) • + + a • (a++a-) . (e^^a+ + e-'S") 

(3.1) 

where A'' denotes the number of sites, af — ^{aj±icr^), and aSR, (pG [0,7r) . This corresponds to 
taking in Eq. (|1.2(l ei = ez = and ea; = acos(0) , ey = — asin(0). The Hamiltonian H[a,4>) (|3.1|l 
for all values of the parameters commutes with the charge operator Q (|2.2|) defined in Sec. 2. 
In order to avoid the appearance of a nonlocal number operator and in complete analogy to the 
known case = (see also Sec. 2), one considers the mixed sector Hamiltonian H{a, (f) H2.13|l 

H{a, 0) = H{a, (/)) ■ Pq + H{-a, (f)) ■ Pi. (3.2) 

Here, the , q€ {0,1}, are the projection operators formerly defined in Eq. (|2.13(l . After 
performing a Jordan- Wigner transformation 



Ck = 





(3.3) 
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(3.4) 



the Hamiltonian H{a, (f>) (|3.2|l reads 

j,k=i L 

where the two NxN matrices B{a,(t>) and B'{a, (f) are given by 
Bj,k{oi,4>) = -5j,fe + 5<5-,+i,fc(l-(5j,jv) + 5^j,fc+i(l-<5fc,Af) 

- f (e*^^,,i4,JV + e-*"*5,, jv4,i) (3.5) 

^j,fc("' |'5j+l,fc(l-(5j,Ar) - \5j,k+i{l~5kM) + f e*'^((5j,i4,JV-(5j,Ar(5fc,i) 

and the fermionic operators Ck and c^, fulfill the anticommutation relations 

{cfe,Cfc'} = {4,4,} = 



The matrices B{a,(j>) and B'{a,(j>) have the properties 

Bt(a,0) = 5(a,0) 

B'^(a» = -B'*{a,(j)) = -B'{a, 



(3.6) 



(3.7) 



and, in particular, are not real except for = or a = which correspond to the previously 
examined case of one (usual) defect pi] . 

In the sequel the Hamiltonian H{a, (f) H3.4|l which is bilinear in the fermionic operators Ck and 
4 l|3.3|l shall be diagonalized to the form 

N 

i?(a,0) =^A,(a,0)a]a^.+£;o(a,0,^) (3.8) 

by means of a Bogoliubov transformation, where the operators at and a\. fulfill the same anti- 
commutation relations as Cfe and c\ (|3.6|l . respectively, and the Aj(a, 0) are real numbers for all 
j = 1,...,A^ due to the hermiticity of the Hamiltonian. EQ{a, cj), N) denotes the ground-state 
energy that depends on the definition of the vacuum. For this purpose we make the ansatz 

N 

ak = ^{gk.jCj + hl^^c]) 

(3.9) 

N 

with fc = 1, . . . , iV and (in general) complex coefficients gk,j , hk.j G C. This generalization compared 
to the transformation used in is enforced by the non- reality of the Hamiltonian H3.4|l . From 

the anticommutation relations for the operators ak and a|, (|3.6|l . one yields conditions on these 
coefficients which are in fact equivalent to the 2N x 2N matrix U defined by 

being unitary. Here, g and h denote NxN matrices with the entries gkj and hkj {k, j = 1, . . . , N), 
respectively. In other words, one has to find a unitary transformation U that diagonalizes the 
Hamiltonian H3.4|) with the additional property 

cue = U* , (3.11) 
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C denoting the 2N x 2N matrix 

f °~ ) , C^ = l2N ■ (3.12) 

\ i-N J 

Note that as a consequence of this property the determinant of U is restricted to the values 
deiU € {1,-1}. 

Necessary conditions for the existence of a transformation as desired are the N equations 

[afe,#(a»] = Afe(a,0)-afc , k = l,...,N, (3.13) 
which, using Eqs. (|3.4|) and H3.9|) . after some calculations result in the equations (fc== 1, . . . , N) 

Afc-fffc = gu-B-hl-B'* 

(3.14) 

Ku-K = gk-B'-hl-B* 

together with their complex conjugated. For convenience, the dependency on the parameters a 
and (j) was dropped and the row vectors gk — {gk,i, ■ ■ ■ , gk,N) and = {hk.i, • . • , hk,N) were used. 
Introducing the notations 

^fe = (^9k,i + hk,i, ■ ■ ■ ,gk,N+ hk,N,gl^i + hl j^, . . . ,glj^ + hl j^^ 

(3.15) 

= 7f (^9k,i-hk,i, ■ ■ . ,gk,N~hk,N,hlj^~glj^, hlj^-g'^j^'j 



and 

r 

^ -B' -B 
with B being hermitian, Eqs. H3.14(l become 

Ak-ik = -i^k ■ B 



B=( Z* ) , (3.16) 



Afe • vj/fe = ■ B . 



One therefore derives the following eigenvalue equations for the squares of the fermion frequencies 
Afe (fc = l,...,iV) 

Kl-$k - ^k-A = $fe.B2 

(3.18) 

Kl-^k = ^k-A = , 
where A = B^ is an — at least — positive semi-definite hermitian 2N x 2 matrix given explicitly 

by 

^^"''^^^U'*(«,0) A*M)) ^'-''^ 

with the two NxN matrices A and A' 

Aj.k{a,(j)) = (5j,fe • (2 - ^^^{Sj^iSk,l+Sj^NSk,N)j - <5j+l,A; • (l-<5j,Ar) - (5j,A_.+i • (l-4,Af) 

+ ae'"l'Sj^iSk,N + ae^'^'I'Sj^NSk,! ~ "2"^'^"' ' (^j,24,Ar-5j,w4,2) 

(3.20) 

Note the following property of Aia, 0) 

CA(a,(A)C = l*(a,^) . (3.21) 



The unitarity condition on U (|3.1U|) leads to the restrictions 

• - • - kk , • - , (3.22) 
or, equivalently, to the unitarity of the matrix U , 

Jl' ^^t'L)-T.U, (3.23) 



i{9-h) {i{g-h)) 

that is buih from the row vectors and ^'fe, A; = 1, . . . , iV. Here, T is the unitary matrix 

T^^i ^.^ ) . (3.24) 



In addition, the fact that U has to comply with Eq. (|3.11() leads to the conditions 

$k-C^^l , ^k-C^^l , fc = l,...,iV , (3.25) 

or, equivalently, 

U-C = U* . (3.26) 

However, in Appendix A it is shown that the latter conditions can be fulfilled by a simple change 
of basis always. Thus it is possible to reconstruct from the knowledge of U the transformation 
U (|XTn|) by C/ = Tt[/ and thereby the Eqs. (|TT7|) . and in fact are necessary and 

sufficient to determine the Bogoliubov transformation H3.9|) . 

In this way the determination of the A^ — up to a sign — has been transformed to the problem 
of computing the eigenvalues of a 2A^ x 2A'^ hermitian matrix and, furthermore, a Bogoliubov 
transformation H3.9|l can be constructed from the knowledge of the eigenvectors. The only thing 
that has been missed so far is the degeneration of the spectrum of A that has to exist if the 
Eqs. I|3.18|l can be solved consistently. This is the purpose of Appendix B, where the origin of this 
property — obviously a symmetry hidden in the equations — is examined. One might think that 
this is not necessary because it follows from the ansatz (|3.9|l , but it is the aim to show that this 
ansatz can always be used and does not lead to equations that cannot be solved. The question 
about the sign of the A^ can be answered as follows. A change of the sign of A^ corresponds to an 
exchange of the corresponding creation and annihilation operators aj. and , which is a special 
Bogoliubov transformation (see also Appendix B), together with a shift in the ground-state energy 
Eq H3.8|1 . In order to let the vacuum (the state annihilated by all annihilation operators) be the 
state with lowest energy, one chooses all Afc to be non-negative. 

Let me stop at this point for a short comment. From Eq. H3.20|l . one recognizes that the matrix 
A' {a, 4>) does not vanish for = as one might have expected. This is of course due to the choice 
of the basis (|3.15|l and it is easy to find a Bogoliubov transformation (|3.9|l that block-diagonalizes 
A{a, 0) (|3.19|l to a form with twice the same real symmetric matrix along the diagonal. In this 
way the problem as usually reduces to the diagonalization of a N x N matrix jJT] and one just 
obtains the whole spectrum doubled. Principally, the block-diagonalization remains possible in 
the general case (fi^O, too, but the transformation may become arbitrarily complicated. However, 
in Appendix B it will be argued that the doubling of the spectrum is a general feature — as it 
must be to fulfill Eqs. (|3.18|) — not limited to the case (f> = 0. 

The eigenvalue equations H3.18|l are now solved by means of an ansatz. This is done the same 
way as in [51 ^2 El by 



(3.27) 



with j = 1, . . . , N, where the eigenvectors $p (and likewise the eigenvalues Ap) are from now on 
indexed by the variable p which will (in general) not be integer-valued. One observes that all but 
six of the 2^^^ linear homogeneous equations summarized in 

A{a, (j)) (3.28) 
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are solved by the relation 

A^=(2sin(|))' . (3.29) 

The remaining six equations then determine the possible values of p through the requirement that 
there has to exist a solution for the variables Vj, fj, Sj, and Sj in the ansatz 13.27|l . This is 
obviously equivalent to demand the coefhcient determinant to vanish. The computation of this 
determinant yields a — highly non-linear — equation for the allowed values of p and therewith 
for the eigenvalues Ap . It is known that in this way one obtains all those eigenvalues Ap whose 
squares are smaller than four 5 . However, this is sufficient to determine the finite-size scaling 
limit, because in that case only low-energy states enter the calculation. 

Of course, one is not able to solve the equation that determines p for any size N of the chain. 
One is rather interested in the thermodynamic limit iV — > cxd. This is performed by expanding the 
determinant in inverse powers of N , where one uses ^2 El 

N -p^Ep + OiN-^) . (3.30) 

The coefficient of the first non- vanishing power in N^^ leads to an equation for Sp that determines 
the scaled eigenvalues in the continuum limit by 

lim [N ■ Ap) = Sp . (3.31) 

The above-mentioned algebraic calculations were performed on the computer. Let me only 
state the result. The final equation for ep{a, (p) reads 

((1 -I- K^) ■ cos(£p) -I- 2k)^ = , K^a cos{(f)) , (3.32) 

and one recognizes that this is just the square of the equation that one obtains for one 'ordinary' 
defect of strength K = acos(0) ^2^]. Therefore the finite-size scaling limit is the same for both 
systems, although they generally differ on finite chains. Thus the operator content that only 
depends on k can be read off from Eqs. H2.20|l and H2.21|l in Sec. 2 directly. 

Therefore, exactly as for an 'ordinary' defect 12, 2], the scaled (relative to the ground-state 

energy Eq{1, 0, N) for periodic BCs) Hamiltonian H{a, (f>) 

H{a,<t>)= lim —■{H{a,cp)~E^{\,0,N) + Ai{a,4>)) (3.33) 

is given by a generator Lo(¥'(a cos((/)))) of a shifted tt(l) KMA 12.18|l . From the positive solutions 
of Eq. H3.32|l one obtains the scaled values for the fermion frequencies (I3.31|) which have to be 
inserted into Eq. (|3.8() . This results in 

#(a,0) = H{n) 

= [^ + \- '^(^) ) +[^+\ + '^('*) ) ^I^'^ + ^^^^ (3-34) 

= LQ{ip{K)) , 

where the shift </?(«) has the same form as for an 'ordinary' defect | 11[I12| (cf. Eq. (|2.21|l ') 

(P{k) — arctan(K) , (3.35) 

4 TT 

with K = acos{(j>) . Here, Eq. H2.15|l with A2(l) = 7rc/6 and the expansion of the ground-state 
energy Eo{a, 0, N) H3.8|l of the Hamiltonian H{a, 0) (|3.4I) 

- Eo{a, ^,N)^Ao-N + Ai(a, </.) + ^fi^ . + o{N-^) (3.36) 
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were used. The 'effective' central charge c(a,(j)) that also depends on K — acos{(j)) only, is deter- 
mined by the shift ip H3.35|) through the relation 

~ 24 ' ^'^■•^^> 

and the fermionic operators ak and bk in Eq. H3.34|l descend from the ak in Eq. 1)3.8(1 by means of 
a suitable re-numeration of the indices. 

The only thing involved in the resulting expressions that does not depend on K = acos{(j)) alone 
is the (non- universal) surface energy Ai{a,4)). Since the ground-state energy does not obtain 
from the analysis performed, no analytic expression is known for this term. However, numerical 
investigations show that it differs from the Ai{acos{(f>)) of Eq. H2.16|l by an amount of i at most 
(for (f)= ^ and a — > oo), thus it resembles the function for an 'ordinary' defect. In analogy to the 
case of 'ordinary' defects, a generalization of this investigation to an arbitrary number of isolated 
generalized defects of this form at commensurate distances ^2 ^| or to a generalized extended 
defect |15j obviously can be performed the same way as for 'ordinary' defects. 



4 Finite-size non-integrable generalized defects 

In this section another two-parameter family of Hamiltonians 7i(/3, V') defined by 

H{/3,^) - -i [^1 + ^J^+i] +< + ■ (cos(V')aJ= - sin(^)aj) j , 



(4.1) 



will be considered. Here, /3eR— {0} and ^€(0, tt), whereby one guarantees that f3sm{'ip) is 
different from zero because /3sin(?/') = recovers finite-size integrable defects examined in the 
preceding section. As has been discussed in Sec. 2, the Hamiltonian 1(4. 1() is not finite-size integrable 
under this circumstances. This is caused by the occurrence of the term (J%crf that not only 
breaks the global Z2 invariance of the Ising chain Hamiltonian (see Eq. H2.2|l 'l. but also leads to 
the appearance of a three-fermion coupling if a Jordan- Wigner transformation as in Eq. 1(3.3(1 is 
performed on H ((4.1(1 . 

Since for this reason the methods used in the previous section can not be applied, the Hamil- 
tonian ((4.1() has been investigated by numerical calculations only. This was done by diagonalizing 
H for up to sixteen sites using an algorithm by Lanczos 25'. The scaled energy gaps Fk{f3, ^p, N), 
defined in this case relative to free BCs 

— N 

FkiP, ^,N) = - {Ek{f3, V', N) + A,i(3, ^) - Eo{0, 0, N) ~ Ai(0, 0)) (4.2) 

TT 

with Ai (0, 0) = — 5 1(2.16(1 . have to be extrapolated to an infinite system 

Tk{P,^)= lim FkiP,^,N) . (4.3) 

Here, Ai{P,tp) denotes the surface term that occurs in the expansion of the lowest eigenvalue 
Eo{l3,^,N) (cf. Eq. (233)) of the Hamiltonian 

-EoW,^j,N) = Ao-N + AMi^) + ^^^^^-N-'+o{N-') , (4.4) 

and Efi.{(3,ip, N) is the fc*'' eigenvalue of the iV-site Hamiltonian T-C{f3,tp) ((4.1() . The partition 
function J-'{P, ip, z) for the finite-size scaling spectrum is defined by 

m^,z)= hm Y^zF^iP^^^^) , (4.5) 
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where the sum extends over the whole spectrum of Ti (|4.1|l . The use of Eo{0, 0, N) as the reference 
energy in Eq. (|4.2|l wiU be justified by the final expression for the partition function ^{(3, ip, z) l|4.5|l 
conjectured on the basis of numerical results. 

There is one disadvantage in the use of Eqs. (|4.2|l and (|4.3|l for the actual numerical compu- 
tations, namely that the analytic form of the surface energy -0) (|4.4|l is not known. There- 
fore, one alternatively extrapolates the scaled energy gaps referred to the relative ground-state 
Eo{f3, tp, N) and besides that determines the lowest scaled energy gap J-a{l3, -0) 

— N 

.Fo(/3,0)= lim - ■ {Eo{(3,ij,N) - Eo{0,0,N) + Ai{f3,^) - Ai{0,0)) (4.6) 

from the finite-size corrections of the relative ground-state energy i?o(/3, '0> H4.4|l . Knowing 
the 'effective central charge' c(j3,tp) 14.4|l is equivalent to the knowledge of !Fo{f3,il;), because by 
Eqs. iPT5|| (with yl2(0) = f|) and (HSJ the relation (c=i) 

ToiM) = 

_ l-2c{p.,-4i) ^ • 

48 

holds. This was used to numerically compute J^o{f3,ip) for various choices of the parameters. For 
this purpose the known value of = 7 was inserted into Eq. H4.4(l and the unknown surface term 
Ai {f3, Tp) H4.4|l was eliminated by combining the equation for two subsequent values of N. The value 
of J^oiP,^) then could be determined approximatively by extrapolating to infinite size N 00, 
wherefore the algorithm developed by Bulirsch and Stoer (5Hlj see also ^3) was employed. 

The numerical results are in reasonable agreement with the following conjecture for the operator 
content expressed through the character functions [20] xa{z) (|2.5|) of the unitary irreps of the c— ^ 
VA H2.4|l corresponding to a highest weight A , Ag{0,j^,^}. For /3>0 one has 

r xo{z) if V^e(0,f) 

^(AV',^) = I XJ,iz) if = f (4.8) 

wherefrom the case < can be deduced using the relation 

T{P,i^,z)^Ti~p,7r-ij,z) . (4.9) 

This equality follows directly from the fact that for arbitrary N the Hamiltonian H4.1|l possesses 
the property 

Q • H(/3, 0) • Q = W(-/3, TT - V) (4.10) 

implying the equality of the finite-size spectra of H(/3, ip) and H{~(3, n — ip). 

In Table 1 the obtained estimates of the lowest scaled energy gap J^oiP, 0) (|4.7|l together with 
their — rather subjective — errors for /3=1 and different values of tp are presented. Table 2 
contains the same thing for three fixed values of ipi namely 7r/5 , 7r/2 , and 47r/5 , and varying p. 
Of course, only a small fraction of the collected data could be displayed therein. In both tables, 
the values deduced from the conjecture 1)4. 8|l are opposed to the numerical data. A comparison 
shows that the agreement is good on the whole, especially if one stays far enough away from 
/3 = and 0G {0, §, tt}. Near those values of the parameters the extrapolation is difficult (and the 
resulting 'errors' thereby big), what can be understood from the assertion that the dependence on 
the parameters ceases to be continuous there. Obviously one is not able to see such an eff'ect in a 
finite sample, thus the observed behaviour agrees with what one would expect. Anyhow one might 
have awaited possible discontinuities at those points where the symmetry of the Hamiltonian (|4.1|l 
changes, these are all the critical values mentioned above with the exception of V'=f • In this 
case no obvious change of the symmetry occurs, but there might be a 'hidden' symmetry that 
one cannot find by inspection. In Table 3 and Table 4, numerical results for the surface energy 
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Table 1: Opposition of numerical estimates and by Eq. H4.8|) conjectured values of the lowest scaled 
energy gap Tq{(3,iP) H4.6|I for j3 — l and different choices of the parameter V'- 





numerical estimate 


conjectured value 


O.Itt 


0.0005(3) 


0.000000 


0.277 


0.0001(2) 


0.000000 


O.Stt 


0.0002(3) 


0.000000 


0.477 


0.001(2) 


0.000000 


0.577 


0.06250(3) 


0.062500 


0.677 


0.6(1) 


0.500000 


0.777 


0.498(1) 


0.500000 


0.877 


0.4995(3) 


0.500000 


0.977 


0.5(1) 


0.500000 



Table 2: The same as in Table 1, but for three fixed values of the parameter ip and varying (3. 



p 




numerical estimate 


conjectured value 


0.5 


0.277 


0.0002(3) 


0.000000 


2.0 


0.277 


0.00003(2) 


0.000000 


5.0 


0.277 


0.00001(2) 


0.000000 


20.0 


0.277 


0.00000(1) 


0.000000 


100.0 


0.277 


0.00000(1) 


0.000000 


0.5 


0.577 


0.06250(1) 


0.062500 


2.0 


0.577 


0.06251(2) 


0.062500 


5.0 


0.577 


0.06250(1) 


0.062500 


20.0 


0.577 


0.06250(1) 


0.062500 


100.0 


0.577 


0.0625(1) 


0.062500 


0.5 


0.877 


0.5(1) 


0.500000 


2.0 


0.877 


0.5005(3) 


0.500000 


5.0 


0.877 


0.50001(2) 


0.500000 


20.0 


0.877 


0.50000(1) 


0.500000 


100.0 


0.877 


0.50003(2) 


0.500000 



Ai{(3,4') H4.4|l for various values of the parameters /3 and ip are given. In this case no analytic 
expression is known and thus there is nothing to compare with. 

Another support for the conjecture 1)4. 8|l results from the counting of the degeneracies of the 
excited levels ('conformal towers'). These could in all cases be seen up to the forth or fifth excitation 
of the lowest level. First, there are no indications of the existence of a second primary field in 
the operator content, that is up to the (at least) forth excited level of the relative ground-state 
nothing but the expected levels — the conformal tower of the ground-state — has been detected. 
On the other hand, the degeneracies of the conformal towers for the three possible primaries at 
c= ^ differ (at the latest) in the third excitation, therefore the observed degeneracies that complied 
with Eq. (|4.8|l in all examined cases, strongly support the conjecture. 

5 Conclusions 

It was the aim of this article to examine the finite-size scaling properties of the quantum Ising chain 
with one 'generalized' defect H1.2|) . which has been introduced by an arbitrary transformation in 
the algebra of observables in the coupling term at one particular site of the chain. The finite- 
size scaling spectrum for the Ising chain equipped with the most universal defect that maintains 
its global 7i2 symmetry was obtained analytically and proven to be the same as in the case of 
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Table 3: Numerical estimates of the surface energy Ai(/3,-0) (|4.4|) for /3 = 1 and different values of 





numerical estimate 


O.Itt 


0.04377(2) 


0.27r 


0.07821(1) 


O.Stt 


0.09682(3) 


0.47r 


0.10143(3) 


0.57r 


0.097454(3) 


O.evr 


0.101(1) 


O.Tvr 


0.095(2) 


O.Stt 


0.080(1) 


0.97r 


0.0434(1) 



Table 4: The same as in Table 3, but for three values of the parameter and varying f3. 







numerical estimate 


0.5 


0.27r 


-0.09663(3) 


2.0 


0.27r 


0.5105577(2) 


5.0 


0.27r 


1.9502580(1) 


20.0 


0.27r 


9.4111164(1) 


100.0 


0.27r 


49.3993176(2) 


0.5 


0.57r 


-0.0798648(2) 


2.0 


0.57r 


0.522629(3) 


5.0 


0.57r 


1.9447786(1) 


20.0 


0.571- 


9.387026(2) 


100.0 


0.57r 


49.368324(1) 


0.5 


O.Stt 


-0.0958(3) 


2.0 


O.Stt 


0.5105(3) 


5.0 


O.Stt 


1.95026(1) 


20.0 


O.Stt 


9.4111165(3) 


100.0 


O.Stt 


49.399318(1) 



one suitably chosen 'ordinary' defect and thus could be expressed in terms of irreps of a shifted 
w(l) Kac-Moody algebra. The result is (see Eqs. (|3.32l) - (|3.35|) ) that for those generalized defects 
that do not break the global Z2 symmetry of the quantum Ising chain the partition function in 
the continuum limit is in fact identical to that for one ordinary defect. The generalization to 
an arbitrary number of defects of this form at commensurate distances obviously is completely 
analogous to the case of 'ordinary' defects ^^Ej- Recently, a similar observation |2HI has been 
made concerning the Ising quantum chain with three spin couplings |28| and the Ashkin- Teller 
model equipped with special defects |29) that break the symmetries in a similar way in both 
systems. This suggests that there is some sort of 'universality' for quantum chains with defects. 
This means that the partition function of the system will not depend on the detailed form of the 
defect, but only on the symmetries of the chain and the symmetry breaking caused by the defect. 

In the remainder of this article another two-parameter family of generalized defects given 
by cr^+i = /? • {cos{ip)af — sm{^)af) with /3gR— {0} and ■0G(O,7r), that destroy the finite-size 
integrability property of the quantum Ising chain has been studied by means of numerical finite-size 
calculations for chains with up to sixteen sites. They suggest a qualitatively different behaviour of 
the finite-size scaling spectrum, which appears to depend non-continuously on the parameter ^ of 
the defect. In this case a conjecture for the operator content involving only one primary field of a 
c = i VA for each choice of the parameters was given and compared to the results obtained from 
the numerical diagonalization of the chain. For some examples of the parameters, this comparison 
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is presented in the four tables, where in addition numerical estimates for the non-universal surface 
energy (cf. Eq. H4.4|l ) are given. For ■!/) = 7r/2, the resulting finite-size scaling spectrum is the same 
as for mixed BCs, whereas for iI}^-k/2 one obtains the finite-size scaling spectrum of one of the 
two sectors of fixed BCs depending on the sign of /3 cos('0) only. Therefore, the highest weights A 
of the c=\ Virasoro irreps are observed to be A = for /3cos(^) (strictly) positive, A = l/16 in 
the case cos('0) =0 , and finally A = 1/2 for negative values of /3cos('0) • 

In this article, not all possible generalized defects were investigated. In principle, one has four 
real parameters ei, e^;, ey, and Cz corresponding to a defect of the form (|1.2(1 for the quantum Ising 
Hamiltonian H2.1|l . so the two types of defects studied in Sec. 3 and Sec. 4 correspond to ei = ez = 
and arbitrary Cx and Cy resp. ei=ey = Q and arbitrary values of Cx and . Thus the question, if 
the finite-size scaling spectra for the chain with a localized defect different from those considered 
herein will show the same behaviour as observed or not, remains unanswered. One might guess 
that this is true, since a localized defect should not affect the conformal properties except that the 
two c=^ VAs describing the translational invariant case become coupled and thereby only their 
direct sum appears in the spectrum generating algebra. From this one would expect that for every 
choice of a localized defect the spectrum generating algebra is at least a (shifted) u{l) Kac-Moody 
algebra. Furthermore the situations of several isolated, but possibly different, generalized defects 
or even extended generalized defects have not been looked upon. Clearing up these points is left 
to future investigations. 
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Appendix A 

In this place a simple, but maybe non-obvious corollary of linear algebra shall be proven. It 
has been used in Sec. 3 to show that one can always find a Bogoliubov transformation (|3.9|l 
diagonalizing the Hamiltonian H3.4I) to the form (|3.8|l and which beyond that is compatible with 
the conditions 13.10|l and 13.11|l . This was necessary to guarantee that the new operators defined 
by Eq. (|3.9(l respect the fermionic algebra H3.6|l . For this purpose one needs only Eq. H3.21|l as 
supposition. This enables us to state the proposition in a more general way. 
Let M be a hermitian 2N x 2N matrix with the property 

CMC^M* , 

where C denotes the matrix (|3.12() 

„ / Oat Iat \ r'Z 1 

\ i-N J 



Then there exists a set B of 2N (row) vectors B = 



{8/c, k—\,..., 2N} fulfilling the three relations 

(A.l) 

Ofe , AfeGR (A.2) 

(A.3) 



= 5j,k 
= Afe. 
= ©fe 



for aUj,fce {l,2,...,2Af}. 

Since M is hermitian, there exists an orthonormal set S'-^-' = {©i°\ k—\,... ^2N} of eigen- 
vectors of M. This means that B'*'^ automatically satisfies the first two relations (A.l) and (A.2) 
of the proposition. The proof is based on an iterative prescription how to built from this set S'"' 
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another set B of orthonormal eigenvectors of M that satisfies the additional third condition in the 
proposition. 

Before commencing with that, consider an arbitrary normahzed eigenvector of M belonging 
to an eigenvalue A, that is 

e-M = A- e , e-et = i , AeR. 

It follows 

(e*c) • M = e* • (CMC) ■ c 

= Q* ■M*-C 

= A-(e*c) , 

i&c) ■ (e*c)t = 1 , 

that is, (0*C) is normalized eigenvector of M to the same eigenvalue A. 

Now, start with 6^ . Then there are two possibilities: either the two vectors (9^ • C) and 
^ that are both eigenvectors of M to the eigenvalue Ai linear depend on each other or not. In 

the first case you have 

0(0)* . c = e'^f . eW 
with 7 e [0, tt). The vector 0i defined by 

has the properties 

@\-c = Qi , ei-el = i . 

Take now as the new set B^^) = {Oi, e^°\ k = 2,..., 2N}, which of course still fulfills (A.l) and 
(A. 2). In the second case, define two new vectors 9i and O2 by 



61 = ^=^fe^^el"' + e-^el")* 



V2(l+d) 

where 



■C 

c) 



with d G [0, 1) and 6 € [0, tt). Now choose 2N—2 new vectors ©[.^^ to be orthonormal eigenvectors 
of M perpendicular to the subspace spanned by Oi and O2 such that the two sets B^^^ and 
BW, = {91,82, e^'\ k = 3,...,2N}, span the same space (one can obviously restrict this 
orthogonalization procedure to the set of vectors Q'f ^ that span the eigenspace of M belonging to 
the eigenvalue Ai). This ensures that the new set B^^^ fulfills the relations (A.l) and (A. 2) and, 
in addition, 

e*-c = ej , ie{i,2} . 

By iteration of this process (that means performing successive changes of basis where in each 
step one or two additional vectors comply with Eq. (A. 3)), one finally — after at most 2Ar steps 
— ends up with a set B that fulfills the proposition. 
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Appendix B 



This Appendix contains another proof that has been omitted in Sec. 3 in order to keep the 
representation easier to survey. It shall be shown that all eigenvalues of the matrix A{a,(f>) H3.19|l 
are at least doubly degenerate (in fact the degree of degeneration is always even). Even more, 
if U (|3.23|) is any unitary transformation compatible with the condition H3.26|l that diagonalizes 
A{a, (j>) for some fixed values of the parameters, then the diagonal matrix Ad = U^A{a, (j>)U has 
the form 

A2 On 



where A denotes the diagonal NxN matrix with elements A^ , k = l, . . . , N. 

The proof of this fact will be organized as follows. First, we will take into account that the 
Bogoliubov transformation H3.9|l is of course not uniquely determined by the claimed properties. 
From this consideration, it will be possible to find two different transformations that yield the 
same diagonal matrix Ad and the conjectured form of Ad follows directly. As will be seen in 
what follows, the real cause of the degeneration lies in the fact that by exchanging the role of 
creation and annihilation operators one changes the sign of the Ak- 

Assume that one succeeded to find a certain unitary transformation U that diagonalizes A{a,(t>) 
for definite values of a and (j) and in addition complies with the condition (|3.26|l . It was proved in 
Appendix A that this is always possible. Then you have 



ADia,^)^U^Aia,q^)U = i ^ , (B.2) 




where now A and A' are diagonal NxN matrices with elements A^ and AJ., respectively 
(A: = 1, . . . , N). Now, suppose you perform a Bogoliubov transformation H3.9(l involving only cre- 
ation and annihilation operators aj and aj with one fixed index j, that is 

(ij = guj + h*aj 

(B.3) 

a] = haj +g*al , 



where the 2x2 matrix u defined by 



9 h* 
h g* 



is unitary (note that this implies that either g or /i is equal to zero and that the determinant 
det(u) G {1,-1}). This transformation will not alter the anticommutation relations, thus the 
operators and dj are fermionic operators as well. On the other hand, the Hamiltonian H (|3.8|l 
given by 

N 

H^V] Akalak + Eo (B.4) 



k=l 



is transformed into 

N 



H = '^{l-Sj^k)Akalak + det{u)Ajd^jdj + Ajhh* + Eq . (B.5) 



fc=i 



This means that if det(u) = l the Hamiltonian is not changed (since this implies h — 0), whereas in 
the case det(u) = — 1 the sign of A^ is altered and the ground-state energy is shifted by the amount 
Ajhh* = Aj. Of course one can perform a transformation as such independently for every value of 

k, k = l,...,N. 

It follows that the equations deduced from the necessary conditions H3.13|l 

[ak,H]^Akak (B.6) 
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are — up to the signs of the — identical, independent of any performed transformation of the 
above-mentioned type. The shift in the constant term Eq (see Eqs. (jB.4|) and (jB.5|) ') drops out 
because of the use of the commutator in Eq. (|B.6p . However, in the final eigenvalue equations H3.18|l 
only the squares of the enter the calculation, therefore these are not altered at all. This leads to 
the consequence that the diagonal matrix is invariant under any of these transformations (|B.3|) , 
where of course one has to take into account the change of basis performed in the derivation of the 
eigenvalue equation. In other words, the whole problem has a J7(l)®^-invariance corresponding to 
an arbitrary phase factor for each pair of creation and annihilation operators. But the eigenvalue 
equations (|3.18|) beyond this are not altered by an exchange of creation and annihilation operators, 
therefore they enjoy an (J7(l)(X)sZ2) -invariance, where iS>s denotes a semi-direct product with 
the normal divisor (invariant subgroup) on the left side. 

Now, take as a special transformation of this form the matrix C defined by 

C=( .^^ ^^'^ ) . 
\ iIn On J 

This acts on the aj^ and ak in the same way for all k, namely through 

Ofc I — > ak = ial 

al I — > al — —iak ■ 
It follows that the transformation V that is related to U — TU (see Eqs. (13.24(1 and H3.23|l ) by 

V = TCU = CTU = CU = CU*C 

with 

^ _ I On \ 7^ — J_ I ■'■^ ■'■^ 1 

^ \ In On J ' ~ V2 y iIn -iIn J ' 

diagonalizes A to the same diagonal matrix Ajj as U does. That means that Ajj has to commute 
with C and, using Eq. I|B.2|I . one therefore obtains the result 

A2 = A'" 

that proves the assertion IIB.l|) . 

References 

[1] E. Ising, Z. Phys. 21 (1925) 613. 

[2] J. B. Kogut, Rev. Mod. Phys. 51 (1979) 659. 

[3] R. J. Baxter, Exactly Solved Models in Statistical Mechanics. Acad. Press, 1982. 

[4] L. Onsager, Phys. Rev. 65 (1944) 117. 

[5] E. Lieb, T. Schultz and D. Mattis, Ann. Phys. 16 (1961) 407. 

[6] A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Nucl. Phys. B 241 (1984) 333. 

[7] D. Friedan, Z. Qiu and S. Shenker, Phys. Rev. Lett. 52 (1984) 1575. 

[8] J. L. Cardy, Nucl. Phys. B 240 (1984) 514, 
Nucl. Phys. B 275 (1986) 200. 

[9] L. Turban, J. Phys. A: Math. Gen. 18 (1985) L325. 



18 



[10] M. Henkel and A. Patkos, Nucl. Phys. B 285 [FS 19] (1987) 29, 
J. Phys. A: Math. Gen. 20 (1987) 2199, 
J. Phys. A: Math. Gen. 21 (1988) L231. 

11] M. Hcnkcl, A. Patkos and M. Schlottniann, Nucl. Phys. B 314 (1989) 609. 

12] M. Baake, P. Chaselon and M. Schlottmann, Nucl. Phys. B 314 (1989) 625. 

13] M. Schlottmann, Diploma thesis, Bonn-IR-88-41 (1988). 

14] M. Baake, P. Christe and V. Rittenberg, Nucl. Phys. B 300 [FS22] (1988) 637. 
15] H. Hinrichsen, Preprint Bonn-HE-89-12 (1989). 

16] T. W. Burkhardt and I. Guim, J. Phys. A: Math. Gen. 18 (1985) L33, 
Phys. Rev. B 35 (1987) 1799. 

17] G. V. Gehlen and V. Rittenberg, J. Phys. A: Math. Gen. 19 (1986) L631. 

18] P. Chaselon, J. Phys. A: Math. Gen. 22 (1989) 2495. 

19] M. Baake., J. Math. Phys. 29 (1988) 1753. 

20] A. Rocha-Caridi, in Vertex Operators in Mathematics and Physics, 

ed. J. Lepowski, S. Mandelstam and I. Singer. Springer, New York, 1985. 

21] M. Baake, Ph.D. thesis, Bonn-IR-87-09 (1987). 

22] H. W. Blote, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56 (1986) 742. 

23] I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 

24] P. Goddard and D. Olive, Int. J. Mod. Phys. A 1 (1986) 303. 

25] C. Lanczos, J. Res. Nat. Bur. Stand. 45 (1950) 225. 

26] R. Bulirsch and J. Stoer, Num. Math. 6 (1964) 413. 

27] M. Hcnkcl and G. Schiitz, J. Phys. A: Math. Gen. 21 (1988) 2617. 

28] T. Wittlich, Preprint Bonn-HE-89-06 (1989), 
Diploma thesis, Bonn-IR-89-34 (1989). 

[29] U. Grimm, Diploma thesis, Bonn-IR-88-30 (1988). 



19 



